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Abstract 

In the Standard Model, b quarks produced in e'’'e“ annihilation at the TP peak have a large average 
longitudinal polarization of —0.94. Some fraction of this polarization is expected to be transferred 
to b-flavored baryons during hadronization. The average longitudinal polarization of weakly decaying 
b baryons, is measured in approximately 4.3 million hadronic TP decays collected with the 

OPAL detector between 1990 and 1995 at LEP. Those b baryons that decay semileptonically and 
produce a A baryon are identified through the correlation of the baryon number of the A and the 
electric charge of the lepton. In this semileptonic decay, the ratio of the neutrino energy to the lepton 
energy is a sensitive polarization observable. The neutrino energy is estimated using missing energy 
measurements. From a fit to the distribution of this ratio, the value 

= -O.bOAilj ± 0.09 

is obtained, where the first error is statistical and the second systematic. 
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1 Introduction 


According to the Standard Model, the process e+e —> > bb gives rise to b quarks that are 

longitudinally polarized ii with a large average value of (P^) = —0.94 for a weak mixing angle of 
s\v? 6w = 0.23. This polarization varies by only ±2% over the full range of the production anglej^ 9. 
If subsequent hadronization to a baryon is considered for example, the light u and d quarks form 
a spin-0 system and the spin of the Aj^ should be carried entirely by the b quark. In the heavy-quark 
limit, an important prediction of Heavy Quark Effective Theory (HQET) is that the degrees of freedom 
of the b quark are decoupled from the spin-0 light diquark so that the A® should retain almost 100% 
of this polarization |p with only a slight reduction of about 3% due to hard gluon emission during 
hadronization |^]. This can be contrasted with spin-0 pseudoscalar B mesons where the polarization 
information of the b quark is lost or vector B* meson where the polarization information is not 
expected to be observable Q . In b baryon production, it is also possible for the b quark to combine 
with a spin-1 uu, ud, or dd diquark system to form the higher mass baryonic states Sb and that 
are expected to decay strongly to A^tt leading to substantial reduction in the polarization averaged 
over the weakly decaying b baryons |^, Measuring the Ab polarization^ therefore provides a test 
of HQET and information about heavy baryon hadronization and nonperturbative corrections to spin 
transfer in fragmentation. 


The sign of measured polarization also gives information on the chirality of the b quark coupling 
to the weak charged current [^, particularly interesting because of recent analyses testing for the 
presence of sizeable right-handed components in the b —> c charged current coupling |^. However, in 
this case, a polarization values of specific b baryon states would have to be assumed to extract limits 
on the right-handed component of the charged current and hence the limits would be strongly model 
dependent. 


In the weak semileptonic decay Ab —>■ Xc^ t'^X, both the charged lepton and neutrino energy spec¬ 
trum are sensitive to {P^^) ill- Assuming that b decays proceed via the usual left-handed current 
(i.e., {V — A) coupling), in the rest frame of the Ab, charged leptons tend to be emitted antipar¬ 
allel to the spin of the Ab and the vis parallel to the Ab spin. In the laboratory frame, polarization 
then implies a harder lepton energy spectrum and softer neutrino energy spectrum compared to the 
unpolarized case. However, there are substantial uncertainties in the exact shape of these spectra due 
to uncertainties in fragmentation, the ratio of the quark masses ruc/mb, and QCD corrections 
such that the extraction of the average polarization from the spectra alone is problematic. These 
systematic effects partially cancel in the ratio of the average of the lepton energy to the average of 
the neutrino energy or in the ratios of higher moments of the energy spectra 12]. 


This paper describes a measurement of the average longitudinal polarization of weakly decaying 
b baryons, using about 4.3 million multihadronic decays collected by the OPAL detector 

from 1990 to 1995. To ensure a data sample with reasonably large statistics, events containing decay 
chains of the form^ Ab —> followed by A^ —> AX are selected without reconstructing the 

intermediate charm state. The correlation of a A with a negatively-charged lepton (£ = e or /r) or a A 
with a positively-charged lepton can indicate the presence of a semileptonic b-baryon decay, as used 

^In the OPAL coordinate system the x axis points towards the center of the LEP ring, the y axis points upwards and 
the 2 axis points in the direction of the electron beam. The polar angle 9, the azimuthal angle (j) and the radius r denote 
the usual spherical coordinates. 

^The symbol Ab will be used in this paper to refer to all b baryons that decay weakly, including for example the 
Ab and the b quark in combination with a spin-1 us or ds diquark forming a Hb. The symbol A® is used to denote the 
particular ground state b baryon with quark content (bud). 

^Charge conjugate processes are implied throughout this paper. 
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previously by OPAL [1^, ]A|. The other charge combinations {A-i~ and are used to characterize 

the background. The neutrino energy was estimated using missing energy measurements in the 
hemisphere containing the K-i~ pair. To extract the polarization of the b baryon, a fit was then made 
to the distribution of the ratio Ey/En formed event by event. Fitting to the shape of this distribution 
should lead to an estimate with lower statistical variance than the use of the average sample energies 
{Ey)/{Efi) as employed in the previously published measurement of |1T|. 


2 The OPAL detector 


A detailed description of the OPAL detector can be found elsewhere p^ ]. The almost complete her- 
meticity of the OPAL detector allows the missing energy in an event to be determined effectively. The 
central tracking detector consists of a two-layer silicon microstrip detector with polar angle coverage 
I cos0| < 0.8 immediately surrounding the beam-pipe, followed by a high-precision vertex drift cham¬ 
ber, a large-volume jet chamber that records the momentum and energy loss of charged particles over 
98% of the solid angle, and chambers which measure the z-coordinate, all in a uniform 0.435 T axial 
magnetic field. Charged particles can be identified by their specific ionization energy loss, dE/dx, in 
the jet chamber. Further information on the performance of the tracking and dE/dx measurements 


can be found in Ref. |n]. A lead-glass electromagnetic calorimeter is located outside the magnet coil, 
which, in combination with the forward calorimeter, gamma catcher and silicon-tungsten luminome- 
ter |l^, complete the geometrical acceptance down to 24 mrad from the beam direction. The magnet 
return yoke is instrumented with streamer tubes for hadron calorimetry and is surrounded by several 
layers of muon chambers. 


3 Monte Carlo simulations 


To model the signal, samples of simulated data were generated of TA hadronic decays that included 
Aj^ —> A^.^“p£X followed by A^ ^ A7r'’“X and A —> pvr using a modified JETSET decay routine 
incorporating polarized Aj^ decay with QCD corrections as described in Ref. [10|. Samples of 40 000 
events each were generated at twelve different values of {P^^). These signal samples and the heavy 
flavor Monte Carlo event samples described below were generated using the latest world average 


branching ratios for Ac decays []R|. Additional samples including primary production and decay of 
Sb or Hb baryons were used to assess the effect of the presence of different baryons on the result. 


Monte Carlo simulation samples of inclusive hadronic decays were used to check backgrounds. 
The JETSET 7.4 parton shower Monte Carlo generator with the fragmentation function of Pe¬ 
terson et al. |^T| for heavy quarks was used to generate samples including approximately 4 million 
hadronic Z^ decays and 2 million ^ bb decays (the equivalent of about 9 million hadronic decays). 
As will be discussed later, in the decay b —> when the r decays leptonically, the resulting lepton 

and neutrino background has to be taken into account carefully. The tau polarization in these decays 


was calculated according to Ref. [^| and leptonic decays of these and all other r leptons were handled 
by the TAUOLA package 


All Monte Carlo samples were processed through a full simulation of the OPAL detector \‘M\ and 
the same event analysis chain was applied to the simulated events as to the data. 
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4 Data and event selection 


This analysis uses data collected at center-of-mass energies within ±3 GeV of the TP mass during 
the 1990-1995 LEP running periods. Following detector performance requirements, track quality 
cuts, and a standard hadronic event selection [^], a data sample of 4.3 million multihadronic events 
was selected. Given the importance of having well contained events in the detector for a reliable 
measurement of the missing energy, each event was required to have | cos 0thr| <0.85 where 0thr is the 
polar angle of the thrust axis. 


The selection of b-baryon decays using A-lepton pairs is similar to the one previously used to mea¬ 
sure the average b-baryon lifetime and production rate |H] but uses improved electron identification, 
a different jet definition, different lepton kinematic cuts, and the inclusion of an additional vr"’" to 
exploit the charge correlation of A-vr'*' in the decay chain Ab —> A+ —> A7r'*“X. These all contribute to 
an overall improvement of the signal-to-purity ratio by 45%. Firstly, electrons and muons having large 
momentum and large transverse momentum with respect to their associated jet (including the lepton 
candidate) were identified. A baryons were then reconstructed and combined with 7r'*"’s that were 
consistent with coming from a common vertex formed by the A and lepton candidate. In the signal 
process the Tr"*" most commonly comes from the A^ decay with a decay length from the point of Ab 
decay which is small compared with the typical A-lepton vertex resolution. Requirements were placed 
on the A-lepton combination followed by requirements to improve the missing energy determination 
and its applicability in estimating the neutrino energy in the semileptonic decay of the Ab. 


Electrons were identified using an artificial neural network which is a simplified version of that 
described in Ref. [^], using only six rather than twelve inputs. These are: the momentum and polar 
angle of the track, the energy-momentum ratio E/p, the number of electromagnetic calorimeter blocks 
contributing to the energy measurement, the ionization energy loss and its error. Photon conversion 
candidates were rejected using another neural network algorithm, similar to that described in Ref. [p^ ], 
but using the new electron identification algorithm to identify the conversion partner to the electron 
candidate. Muon candidates were identified by associating track segments in the outer muon detectors 
with tracks extrapolated from the central tracking detectors \‘2/V\. 


Charged tracks and electromagnetic clusters not associated with charged tracks were combined 
into jets using a cone algorithm |^] with a minimum jet energy of 5.0 GeV and cone half angle of 
550 mrad. Electron and muon candidates were required to have momentum greater than 4.0 GeV and 
transverse momentum with respect to the associated jet axis, pt, greater than 0.7 GeV to obtain a 
sample enriched in b-quark events. 

Candidate A baryons decaying via A —> pvr” were identified as described in Ref. [1l4| by considering 
oppositely charged tracks assigned to the same jet as the lepton and satisfying AE/dx criteria and 
various kinematic and geometric requirements. To reduce Kg —> 7r'*“7r“ contamination, the invariant 
mass, measured assuming a pion mass for both particles, was required to be more than one standard 
deviation in mass resolution away from the Kg mass: less than 0.491 GeV or greater than 0.503 GeV. 
To decrease the contribution of A candidates coming from fragmentation processes, the A candidate 
was required to have a minimum momentum of 4 GeV. The invariant mass of the reconstructed A 
was required to be between 1.1078 and 1.1234 GeV. The fraction of A baryons wrongly identified as 
A is negligible since for these A momenta the proton candidate can always be taken as the higher 
momentum track. The efficiency of this selection varies from 23% to 10% for low to high momentum 
A candidates. 
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Lepton and A candidates were then correlated using the following criteria. The reconstructed 
momentum vector of the A candidate was required to be closer than 50° to the momentum vector of 
the lepton candidate. The invariant mass of the A-lepton system was demanded to be greater than 
2.2 GeV to reduce the background from A+’s that decay semileptonically: A+ —> Af'+t'^X. To reduce 
further random combinatorial background, the magnitude of the momentum vector sum of the A and 
lepton had to be greater than 9 GeV. 

For each accepted A-lepton pair, candidate 7r'''’s were selected among all charged tracks in the jet 
to which the lepton was assigned. Gharged tracks with momentum greater than 0.5 GeV and with 
measured dE/dx having a probability greater than 1% of being consistent with a pion were considered 
for combining in the vertex. The track was required to be positively charged if combined with a A, and 
negatively charged if combined with a A. A decay vertex was formed in the r-tj) plane by extrapolating 
the candidate A momentum vector from its decay vertex to the intersection with the lepton candidate 
track. The tt candidate was then also fitted to the A-lepton decay vertex and the candidate track with 
the largest corresponding vertex probability was chosen. Ki~Tr~^ combinations were retained if the 

probability of the common fit vertex was greater than 1% and the invariant mass of the common 
vertex less than 5.6 GeV. 

Even for properly reconstructed A.I combinations coming from the decay Ab —> the sub¬ 

sequent semileptonic decay of A+ —> A^+t'^ would result in a second neutrino compromising the 
measurement of E^, from the Ab decay using the missing energy in the event. Therefore a small frac¬ 
tion of the combinations were rejected if the tt"*" candidate track selected above had momentum 
greater than 2.0 GeV and was identified as an electron or muon. 

Each event was divided into two hemispheres by a plane perpendicular to the thrust axis. The 
hemisphere containing the combination was defined as the signal hemisphere and the opposite 

hemisphere as the recoil hemisphere. To reject events with incorrectly measured total energy and 
therefore fake missing energy due to mismeasurement of momenta, signal hemispheres were rejected if 
they contained high-momentum tracks close to anode planes in the jet chamber, i.e., with measured 
momenta greater than 10 GeV and satisfying \(j) — ^anodel < 1° where (j) is the azimuthal angle of the 
track and (/>anode is the angle of the closest anode wire plane in the jet chamber. In the hemisphere 
opposite to the signal hemisphere, additional missing energy can arise due to semileptonic decays of b 
and c hadrons into electrons and muons. Signal hemispheres were rejected if the opposite hemisphere 
contained an electron or muon satisfying the kinematic cuts described earlier since the presence of 
a lepton can indicate the possible presence of a neutrino with substantial energy in the opposite 
hemisphere from a semileptonic decay. 

The invariant mass distribution of the pvr” combination is shown in Fig. Q for Ai~7r~^ (right-sign) 
and A^"’“7r'’“ (wrong-sign) combinations after all other selection requirements have been applied. A total 
of 912 right-sign and 316 wrong-sign A^^vr'*' combinations are selected with an overall b baryon purity 
of 69% in the right-sign sample. The 596 excess right-sign combinations can be attributed mainly to b- 
baryon decays. The overall efficiency to select true A’s and leptons from b baryon semileptonic decays 
is estimated from Monte Carlo simulated samples to be approximately 5.6%, although a knowledge of 
this efficiency is not necessary for the polarization measurement. The number of observed wrong-sign 
combinations will be used after small corrections as an estimate of the level of the background as well 
as to estimate the shape of the background in the distribution of the ratio Ey/Eg. 
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m(p7r-) (GeV) 


Figure 1: Invariant mass distribution of p7r“ combinations in the right-sign A£“7r^ sample (open 
histogram) superimposed by wrong-sign A^+vr"*" (shaded histogram) combinations. Charge conjugate 
states are implied if not stated otherwise. The signal region is shown by the bracketed range of masses. 
The peak at the A mass in the wrong-sign distribution is indicative of genuine A particles (though not 
necessarily from b-baryon decay, e.g. from fragmentation) being combined with real or fake leptons 
(i.e., hadrons misidentihed as leptons). 

5 Neutrino and lepton energy measurement 


In the semileptonic decay of the Ab, the neutrino energy was estimated by the measured missing 
energy in the signal hemisphere, which can be defined by 

_ phemi _ p _l ZT" _ ZT'hemi 

— -^miss — -^beam > -^corr -^vis ’ 

where F^beam is the beam energy and the visible energy in the signal hemisphere. If we assume 

a two-body decay of the Z^, with each hemisphere considered as one body, then Ficorr is a correction 
term calculated assuming that the two hemisphere momenta are equal: 

Ficorr = (iiFgjg — Afj.g^Qjj)/4F^beam- 

Msig and Mj-ecoii are the measured invariant masses of the signal hemisphere and opposite recoil 
hemisphere, respectively. The term F^corr uses the beam energy constraint to improve the resolution 
for the missing energy in the signal hemisphere by accounting for fluctuations of energy splitting 
between hemispheres and unobserved neutral energy on an event by event basis. 

Calculations of these quantities used the four-momenta of reconstructed charged tracks and of 
clusters in the electromagnetic and hadronic calorimeters not associated with charged tracks. The 
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masses of all charged particles were set to the charged pion mass except for the proton from the 
identihed A and the identihed lepton which are set to their appropriate masses. The invariant masses 
of the calorimeter energy clusters were assumed to be zero. Calorimeter clusters associated with 
charged tracks were also included after the expected calorimeter energy for the associated track is 
subtracted from the cluster energy to reduce double counting |^|. If the energy of a cluster was 
smaller than the expected energy for the associated tracks, the cluster energy is not used. 


The agreement of the quantity between data and Monte Carlo events containing identified 

leptons following a preselection with relaxed p and pt cuts to enhance statistics is shown in Fig. |^a). 
To further test and calibrate the measurement of this important quantity, event samples from both 
data and Monte Carlo simulation were prepared that were enhanced in b quarks using a lifetime tag 
in one hemisphere. Hemispheres were identified as containing b hadrons by reconstructing secondary 
vertices and using the parameters of the vertices in an artificial neural network trained to reject non-b 
background [^]. The most important inputs to the neural network were the decay length in the x-y 
plane, its uncertainty, and the number of tracks in the vertex. The selected sample had a b-quark 
purity of approximately 90%. 


Two control samples were used: one with an electron or muon candidate identified as previously 
described in the hemisphere opposite the b-tagged hemisphere and another with no lepton identified in 
the hemisphere opposite the b-tagged hemisphere. The first sample is enriched in energetic neutrinos 
while the second sample is depleted. From the distributions of the measured and true value of 
in these Monte Carlo samples, small residual neutrino energy shifts were observed. These could be 
approximated well with an additional average linear correction, determined from both control samples, 
subtracted from the reconstructed value of that ranged from 50 MeV for = 0 GeV (i.e., 

no neutrino) to 910 MeV for = 30 GeV. This correction was applied to both the data and Monte 
Carlo samples to eliminate the small bias inherent in the neutrino energy reconstruction method. The 
agreement between data and Monte Carlo simulated events in these control samples can be seen in 
Fig. |^(b) and (c). Independent samples of events with hemispheres failing the b-tag requirement, both 
with and without identified leptons, also show good agreement between data and Monte Carlo. 


The absolute resolution on the missing energy and hence the neutrino energy improves with in¬ 
creasing neutrino energy from a{E^) = 5.0 GeV aX Ey = Q GeV to 2.9 GeV at Ey = 15 GeV. In the 
signal Monte Garlo samples, the average neutrino energy resolution was observed to be 3.5 GeV as 
shown in Fig. ^(d). The average neutrino energy for the unpolarized signal process is about 6.1 GeV 
with an r.m.s. spread of 4.4 GeV after applying only kinematic and geometric requirements but not 
including detector resolution effects. 

For the electron and muon candidates the momentum was used as an estimate of the lepton energy. 
For high energy electrons, the electromagnetic cluster energy associated with the electron may have 
better resolution, but the momentum was used to avoid biases due to overlapping energy deposits of 
nearby particles. Changes in the electron spectrum due to final-state radiation are included in the 
Monte Carlo simulation. The distributions of lepton spectra in the control samples described above 
were also compared between data and Monte Carlo. The shapes agree well and the average values are 
{Ef) = 10.18 ± 0.03 GeV for the data and {E() = 10.11 ± 0.03 GeV for the Monte Carlo simulation 
for the b-tagged sample containing identified leptons. 


The distributions of the measured electron and muon energies and of the reconstructed neutrino 
energies for the selected right-sign and wrong-sign A^+Tr"*" samples in the data are shown in 

Fig. |. 
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Figure 2: (a) Reconstructed missing energy in data (points with error bars) and Monte Carlo simulated 
events in qq events after a preselection with relaxed p and pt requirements; the open histogram shows 
the missing energy for events containing h and c hadrons decaying semileptonically, and the shaded 
histogram is for all other processes. Reconstructed missing energy for tagged bb events in data (points 
with error bars) and Monte Carlo simulated events (open histograms) (b) with an identified lepton in 
the hemisphere opposite the b-tag and (c) with no identified lepton in the hemisphere opposite the 
b-tag. (d) Neutrino energy resolution from reconstructed and true in Monte Carlo signal 

events with —> A'^i~uiN. 


6 Background estimate 


The wrong-sign sample is used as an estimate of the level and shape of the background in the right-sign 
sample. The component of the background due to random combinations of tracks and fakes due to 
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Figure 3: (a) Lepton (electron and muon) and (b) reconstructed neutrino energy spectra for the 
hemispheres containing selected right-sign A£~Tr^ combinations with the spectra for wrong-sign com¬ 
binations overlaid. 


misidentified leptons or A baryons is expected to contribute equally to right- and wrong-sign samples, 
as verified by Monte Carlo simulations. 

A larger source of background is the combination of a A baryon from the fragmentation process 
with a genuine lepton from the semileptonic decay of b or c hadrons. In the framework of string 
fragmentation, there is a strong correlation between, for example, the Ab subsequently decaying as 
Ab —> and a high momentum A as one of the neighbouring particles in the fragmentation 

chain. In this case, the background of a fragmentation Ai~ combination will preferentially enter the 
wrong-sign sample. Table gives the background composition of the right- and wrong-sign sample 
as predicted by Monte Carlo simulation and shows the resulting imbalances. These predictions use 
the “popcorn” model |32] and the indicated uncertainties include modelling effects as described in 
detail in Ref. |M|. Note that even if a true lepton from a semileptonic Ab decay is combined with 


a fragmentation A, the ratio Ey/Ei is still indicative of b-baryon polarization since missing energy 
and hence the neutrino energy is calculated using event quantities. Events of this type are therefore 
treated as signal in the fitting procedure (see below). 


A number of exclusive channels as described below are expected to contribute to the right-sign and 
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Source 

RS 

WS 

Signal; A and £ from b baryon 

629 ± 30 


Fragmentation A plus: 

£ from B meson 

94 ± 9 

77 ± 8 

£ from c baryon 

9 ± 3 

3 ± 2 

£ from b baryon 

5 ± 4 

66 ± 6 

£ from c meson 

32 ± 6 

43 ± 7 

Exclusive Backgrounds: 

32 ± 7 

27 ± 6 

Combinatorials and Fakes: 

111 ± 8 

100 ± 8 

predictions of the composition among the 912 right- 


Table 1: Monte Ci 

and the 316 wrong-sign combinations. Errors are due to Monte Carlo statistics, branching ratio 
uncertainties, and modelling systematic errors (see Ref. 


wrong-sign samples, but the total contributions are predicted to be small compared to the background 
due to combinatorials and fakes and fragmentation A baryons. For example, a b baryon can decay 
semileptonically into a r lepton that subsequently decays into an electron and muon to enter into 
the right-sign sample as background. The r polarization in this case was varied within its range of 
uncertainties [^] allowing the At, polarization to vary between 0 and —1 to estimate the systematic 
error on the Monte Carlo prediction for this rate. The expected contribution to the wrong-sign sample 
for leptons arising from the process A+ —> is kept small by the cut on the A£ invariant mass. The 

uncertainty in this rate can be affected by uncertainties in the polarization of the Ac which is varied 
by reweighting Monte Carlo simulated events as described in Section 9. Other exclusive backgrounds 
include small fractions of Ab —*■ A^D“ followed by D“ —> £~X and B —> YcX£i?, (where Yc denotes 
any c baryon) H]- 


7 Fitting procedure 


The b-baryon polarization was extracted by comparing the reconstructed distribution of Ey/En in the 
data to spectra estimated from fully simulated Monte Carlo events corresponding to various values of 
^pAb) satisfying all the selection criteria. A binned maximum likelihood fit was used to extract 
^pAb) determining which Monte Carlo spectrum gives the best description of the data. To obtain 
the Monte Carlo spectrum for arbitrary values of polynomial fits are made bin by bin as a 

function of {P^'°) to allow interpolation between spectra. The binned likelihood method has several 
advantages over analytical functions used to describe the data distributions and resolutions in an 
unbinned likelihood fit: it takes into account correlations between E^, E^, a{Ey) and a^E^); effects of 
kinematic cuts such as the minimum lepton energy; energy, momentum, and missing energy resolution 
with a full simulation of the detector response; and any variation of the selection efficiency with Ey/En. 

The probability density function (PDF) was estimated by the binned and normalized Monte Carlo 
spectrum with being the value of the PDF in bin i of 30 bins for the predicted signal 

distribution for the average polarization value (P^^). is the normalized Monte Carlo spectrum 

of the wrong-sign backgrounds listed in Table |l] excluding the contribution of a fragmentation A 
combined with a genuine lepton from a b-baryon decay. The expected number of entries in bin i of 
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the right-sign and wrong-sign distributions are and given by: 


re, 


ws 




ws 


= 


(1 - /WS)pSig^pAb^ ^ jWSpWS^pAb^ 


jAb 


WS-nWS^oAb' 


■N, 


RS 


tot ) 


withpWS(pAb) ^ (1-/WS^P| 


back 


+ fZ''PTHPh 
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Here is the estimated fraction of fragmentation A combined with leptons from b baryon decays 
as given in Table |l], and are the total number of right-sign and wrong-sign combinations 

observed in the data, and 


The binned likelihood C is the product of Poisson probabilities for obtaining the numbers of 
events observed in the data in each bin using the Monte Carlo expectation in that bin for a given 
(pZ) for both the wrong-sign and right-sign spectra. A normalization constraint was imposed on 
the total number of observed right-sign combinations. An adjusted likelihood [^] was used to take 
into account statistical fluctuations in the Monte Carlo prediction. The above assumes that both 
the level and shape of the wrong-sign sample can be used to estimate the background in the right- 
sign sample. After removing the signal component in the Monte Carlo wrong-sign combinations, the 
distribution of EyjEf^ is consistent with the distribution of the background in the right-sign Monte 
Carlo combinations. Effects of differences are addressed as systematic errors below. 


Further checks of the fitting procedure were performed using Monte Carlo samples. Fully simulated 
Monte Carlo subsamples of known each corresponding to the same number of combinations 

observed in the data were used as input to the fit and no bias was observed at the precision stud¬ 
ied. A simple Monte Carlo program was used to generate an ensemble of 2000 samples, each with 
the same statistics as selected in the data. For a given value of lepton and neutrino energies 

were randomly sampled from analytical joint energy distributions, smeared according to parameter¬ 
ized energy resolution functions, and kinematic cuts applied. The background and distribution was 
generated assuming zero polarization and the wrong-sign distribution sampled separately, including 
a signal component as expected in the data. In all cases, no significant additional systematic biases 
were observed, and residuals divided by the error on {P^^} indicated a correct evaluation of statistical 
errors by the likelihood fit. Lastly, the complete analysis was applied directly to the approximately 
4 million fully simulated hadronic Monte Carlo events. The fitted value of (P^) — 13 is 

consistent with a value of zero polarization for b baryons as simulated in the Monte Carlo sample. 


8 Fit result 


The event-by-event distribution of EyjEi was formed, as shown in Fig. 0(a), and the right-sign and 
wrong-sign distributions were fitted using the described procedure. The resulting curve of — log C ver¬ 
sus (P^^), offset so the maximum value of log C is zero, is shown in Fig. §(b) indicating a measurement 

(P^) = -0.561°:?°, 

and a value of zero polarization ruled out at the 95% confidence level (CL), considering only the 
statistical errors. The per degree of freedom is 0.91. 
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Figure 4; (a) Bottom: distribution of the reconstructed neutrino energy E^, divided by the electron 
or muon energy Ei for each right-sign Af~7r~^ combination in the data (solid circles) overlaid by the 
distribution for each of the wrong-sign combinations (open circles). The solid and dashed lines show the 
result of the likelihood fit to the right-sign and wrong-sign distribution, respectively. Top: fractional 
difference (points with error bars) compared to the analogous distributions 

for the prediction for {P^^) = 0.0 (dotted line) and {Pi'") = —1.0 (dashed line), (b) Negative log 
likelihood as a function of {P^'^) for the fit. The dashed and dotted lines indicate the likelihood change 
required for one standard deviation errors and 95% CL limits, respectively (statistical errors only). 
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9 Systematic uncertainties and consistency checks 

The sources of systematic error that have been considered are summarized in Table 


Source of Uncertainty 


Ey resolution 

±0.02 

Ey reconstruction 

±0.05 

El scale and shape 

±0.03 

Selection criteria 

±0.02 

Background fraction and shape 

±0.04 

b fragmentation 

±0.03 

Ac polarization 

±0.02 

b —> T 

±0.01 

Fitting method 

±0.03 

Theoretical uncertainty (form factor 

±0.03 

modelling, QCD corrections, ruc/mb) 


Total 

±0.09 


Table 2: Summary of systematic uncertainties in the measurement of 

The overall rate or efficiency for tagging A^vr'*' combinations is not used in the fit but any discrep¬ 
ancy in rate between data and Monte Carlo prediction or distortions as a function of Ei,, E^, or their 
ratio can potentially affect the result. Control samples were used in the following to estimate these 
effects. 

To test the Monte Carlo modelling of reconstructed neutrino energy resolution, events with missing 
energy in the b-tagged Monte Carlo control samples with identified leptons in the range corresponding 
to the missing energy in the Ab ^ Eclui signal were studied. Comparing the resolution a{Ey) between 
the signal sample and the control sample in bins of Ey, differences of up to 12% were observed. 
Propagating this resolution uncertainty into the PDF’s for Ey/E^ resulted in differences of ±0.02 in 

The linear correction to used to reconstruct Ey was varied within its uncertainties in slope 

and intercept resulting in an observed variation of ±0.03 in Events with a lepton opposite 

a hemisphere failing a b-tag were used as an independent data control sample. In this sample, the 
data and Monte Carlo distributions were normalized to the same number of events and the ratio of 
the resulting distributions in bins of was fitted to a first-order polynomial. The fitted 

slope was consistent with zero, but combinations were reweighted according to the one standard 

deviation bounds on a possible slope. It was checked that a fit to a second-order polynomial resulted 
in no significant parabolic coefficient. Similar results were observed with a sample consisting of events 
with a lepton opposite a hemisphere passing the b-tag. Variations ranged up to ±210 MeV in the 
average value of Ey and resulted in a change of ±0.05 in (P^^). This is assigned as the systematic 
error due to uncertainties in the Ey reconstruction and its distribution. 

The absolute scale of the momentum measurement of electrons and muons was determined by com¬ 
paring e“’'e“ and pair and photon conversion data with corresponding Monte Carlo simulation. 

Lepton momenta in the data were rescaled by the observed uncertainty of ±0.4% resulting in a change 
in (Pi^) of ±0.01. Possible mismodelling of the electron identification efficiency as a function of Ee 
was investigated by comparing the relative rate in multihadronic events of cleanly identified photon 
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conversions with Monte Carlo predictions as a function of electron energy. Similar studies |^| were 
done for muons as a function of muon momentum. In both cases, relative variations did not exceed 
3%, and Monte Carlo E(^ spectra were reweighted according to the observed differences. Refitting 
resulted in observed changes of ±0.03 in The uncertainty in the momentum resolution results 

in a negligible effect on the polarization and the total assigned systematic error due to uncertainties 
in the distributions was ±0.03. 


The selection criteria were varied within typical ranges of their resolution as determined by signal 
Monte Carlo samples and the fit repeated. The observed variations of ±0.02 in {P^'°) were assigned 
as a systematic error. 


The effect of uncertainties in the determination of the shape of the wrong-sign background PDF 
and the level of the background was assessed by varying each component by its uncertainty indicated 
in Table |l|. Since the background due to a fragmentation A baryon combined with a lepton from 
b-baryon decay was predicted to be (21 ± 3)% of the wrong-sign sample, the fit was also repeated 
after scaling down the background in the right-sign sample by this factor. Lastly, instead of using the 
wrong-sign background PDF fit to the number of wrong-sign combinations, the distribution of Ey/E^ 
in the data wrong-sign combinations was simply subtracted from the distribution of the right-sign 
combinations and the fit to the signal PDF’s repeated. The observed changes from the central value 
of in these tests were added in quadrature resulting in a systematic error of ±0.04 assigned due 

to uncertainties in the background estimate. 


The fragmentation parameter (x£;)b, defined as the average fraction of the beam energy carried 
by a weakly decaying b hadron, has been estimated to be the same for baryons and mesons at the 1% 
level 111]. Varying the average energy by the measured errors on the average b-hadron energy |^] 
added in quadrature with this 1% uncertainty results in a variation of (P^'’) of ±0.03 due to the 
residual logarithmic dependence on fragmentation of the Ey/E^ distribution 


Polarization of the A+ baryon can affect the acceptance for combinations in the signal chain 
Ab ^ A+e-ue, A+ ^ AX, A —> pvr and also the lepton spectra from background semileptonic decays 
of Ac either produced in the above chain or via direct production in —> cc events. In the first two 
cases, measured Ac decay parameters were used and A^tt'’' Monte Carlo combinations involving Ac 
baryons entering the sample for the background PDF were reweighted (also see Ref. to correspond 
to variations of the Ac polarization between 0 and —1. In the latter case, the Ac polarization was 
varied between 0 and —0.68, the value for c quarks produced in —> cc. The fitted value of (P^^) 
changed by ±0.02, which is assigned as a systematic error due to this effect. 


The measured branching ratio for Br(b —> rP^-X) was varied within its uncertainties and the 
level of transfer of polarization of the tau chosen from an alternative model [^| to reweight Monte 
Carlo background combinations. The effect was small, and a systematic error of ±0.01 due to this 
background source uncertainty was assigned. 


Monte Carlo tests for possible biases in the fitting method have already been described earlier. In 
addition, the binning of the spectra and the PDF’s was varied from 30 bins to 20 and 40 bins. The 
range of the fit was also changed from —2 < Ey/E^ < 4 to — 1 < Ey/E^ < 3 and the fit repeated. 
Adding the observed variations in (P^*’) in quadrature, a systematic error of ±0.03 was assigned due 
to possible biases introduced by the fitting method. 


Many sources of uncertainty in the lepton and neutrino energy spectra partially cancel in the ratio 
Ey/Ef. For the input to the theoretical calculation of these spectra for polarized b baryons including 
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QCD corrections ||lOf| , the ratio mdraY^ was varied in the range 0.22-0.39 and the strong coupling 
constant 05 (mb) in the range 0.20-0.24 


An alternative simple quark model |37] using different 
form factors in the polarized Aj( decay was also substituted. The signal Monte Carlo spectra were 
weighted accordingly. Different masses of b baryons and c baryon decay products (e.g., Hb and He) 
were also used in the calculation of the spectra. In all cases the fit was repeated and from the observed 
variation in a systematic error of ±0.03 was assigned due to theoretical uncertainties in the 

predicted E^/Ei spectra. 


The systematic uncertainties from each source were added in quadrature to obtain an estimated 
total systematic error of ±0.09 on 

Consistency Checks: Further checks were performed to search for other systematic effects. The 
cut on the minimum momentum of the electron or muon was varied between 3.0 and 5.0 GeV, and 
the requirement on the minimum pt changed between 0.6 and 1.0 GeV. These variations change the b 
quark content of the sample; in particular the Ac and charm meson content increases for the smaller pt 
cut. The minimum momentum of the A candidate was varied in the range 3.0 to 5.0 GeV, increasing 
and decreasing the background level of fragmentation A baryons being combined with leptons. In each 
case the analysis was repeated and (Pi^) redetermined. Observed variations were consistent with the 
statistical errors on the uncorrelated fractions of the different samples. 


The fit was performed separately for data collected in 1990-1993 and for 1994-1995, and also 
separately for Ae^vr'*' and A/x^tt'*' combinations. All results were statistically consistent with each 
other and with the central value from the full data set. The observed rate of Al'“7r+ combinations 
were also consistent in these different samples, and the derived value /(b ^ Ab) • Br(Ab —> AluiX.) = 
(2.78 ± 0.16) X 10“^ (statistical error only) agrees with previous OPAL measurements ^ within 
errors. 


A fit was also made to the reconstructed E^/Ei distribution in the inclusive lepton control sample 
from the data. This large statistics sample consisted of events with an identified lepton in the hemi¬ 
sphere opposite to a b-tagged hemisphere. No wrong-sign background subtraction was performed, 
and a longitudinal polarization of {Pl) = —0.043 ± 0.019 (stat.) was measured. Although B mesons 
are expected to be unpolarized (the ground states are unpolarized and for spin-1 states the polar¬ 
ization is undetectable), a fraction ( 10 .ll 3 ®)% |[^ of the b hadrons are predicted to be b baryons 
in this sample, and there will also be a contamination of polarized c baryons that decay weakly. In 
addition, the smaller expected semileptonic branching ratio for b baryons compared to the other b 
hadrons ||3^ , [l^ would result in further reduction of total observed polarization. Although the total 
observed polarization is not necessarily expected to be zero, it should be small, as observed. 

A further cross check was made with a data sample containing candidates for the exclusive decay 
B*^ —> 'D*~£^vi selected from events where a and a lepton of opposite charge were found in the 
same jet |^. A value consistent with zero, {Pl) = 0.071qj 2) found as expected for a sample 
with a b hadron content almost exclusively from B mesons and much less contamination from possibly 
polarized heavy baryons. 


10 Discussion and summary 


The average polarization of b baryons in decays at OPAL has been measured to be: 

(P^) = -0.56l°;20±0.09, 
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where the first error is statistical and the second systematic. This level of polarization is larger 
than, but consistent with, the published measurement by the ALEPH Collaboration of (Pi^) = 
(stat.) (syst.). 

Including systematic errors, this OPAL measurement implies bounds on the longitudinal polariza¬ 
tion of b baryons of —0.13 > (P^^) > —0.87 at 95% CL, therefore disfavoring full observed average 
polarization of —0.94. This is the first measurement to exclude zero polarization of b baryons at 
larger than 95% CL providing direct evidence that the b quark is longitudinally polarized in the decay 
ZO ^ bb. 


A simple model Q can be used to predict the total observed b baryon polarization after depo¬ 
larization of those b baryons proceeding through intermediate states involving strong decays: b — 
^ Ab. The prediction depends on a parameter A that is the relative probability of producing a 
spin S' = 1 diquark as opposed to a S = 0 diquark, and a parameter uji that is the probability that 
the spin 1 diquark has angular momentum component ja = ±1 along the fragmentation axis. For 
(Ji = 0 (alignment suppression) and the default JETSET parameteiQ for A, a value of {P^'") = —0.68 
is predicted. Varying u>i to 0.66 (isotropic diquark spin distribution) and A within bounds suggested 
by measurements |41] gives a range of predictions between —0.54 and —0.88. In this model, the mea¬ 
sured value of {P^^) is consistent with no depolarization during fragmentation. It is also consistent 
with the inclusive A measurement of OPAL |^] that also found no evidence of polarization loss in 
hadronization. 
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